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Abstract 

The master equation of one-dimensional three-species reaction-diffusion 
processes is mapped onto an imaginary-time Schrodinger equation. In many 
cases the Hamiltonian obtained is that of an integrable quantum chain. Within 
this approach we search for all 3-state integrable quantum chains whose spectra 
are known and which are related to diffusive-reactive systems. Two integrable 
models are found to appear naturally in this context: the C/ g 5C/(2)-invariant 
model with external fields and the 3-state UqSU(P/M)-mvariaxit Perk-Schultz 
models with external fields. A nonlocal similarity transformation which brings 
the Hamiltonian governing the chemical processes to the known standard forms 
is described, leading in the case of periodic boundary conditions to a general- 
ization of the Dzialoshinsky-Moriya interaction. 
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1 Introduction 



Since the pioneering work of Smoluchowski in 1917 |TJ], reaction-diffusion-limited 
processes have had a forefront position in nonequilibrium statistical Physics. They 

k 

can be portrayed as bimolecular processes of the type A+B # C+D where molecules 

r 

of species A and B (C and D) or two different states of the same molecule react to 
form C and D (A and B) with a reaction rate k (r). Particularly in the last decade 
a great amount of research has been reported (see @-0 and references therein). 
The interest has been centered mainly on irreversible (r = 0) and vacuum-driven 
chemical reactions, i.e. those for which at least one of the final products is an inert 
state (a precipitate or a non-reacting molecule) denoted by 0. The commonly studied 
reactions are 



1. Diffusion 



A + 
5 + 



+ A 
+ B 



2. Interchange 



3. Death 



A + B # B + A 



A + 
B + 



+ 



4. One- and Two-Species Annihilation pO 



A + A 
B + B 
A + B 



+ 



5. Coagulation [111 



6. Trapping 12 



7. Mutation 



A + A -> A + 
B + B -> B + 



A + B -»■ A 
A + B -> £ 



A + -> £ 
5 + -> A- 



1.1) 



;i-2) 



(1.3) 



;i.4) 



(1.5) 



;i.6) 



(1.7) 
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8. Transmutation 



A + -> + 5 

5 + -> + A (1.8) 



9. Polymerisation [fl3|1 



A + A -> 5 + 

B + B -> A + (1.9) 

The apparent simplicity of the processes depicted above is quite deceiving. From 
the great variety of nonequilibrium problems, reaction-diffusion processes are one of 
the most difficult to tackle. An important step towards circumventing these diffi- 
culties, largely of a mathematical nature, was taken by Glauber in 1963 ]14| and 



subsequently explored by many authors [|15 |. In a pursuit to understand nonequib 



rium systems in terms of the more treatable equilibrium ones, he devised an ingenious 
way of using classical spin systems to study the problem of critical dynamics by means 
of a master equation approach. This opened the possibility of employing results from 
spin chains in the context of nonequilibrium problems. The status quo remained how- 
ever pratically unaltered. This was so since the developments achieved on the theory 
of spin chains (or more generally speaking quantum chains), which ultimately lead 
us into the concept of integrability fl7l , were made quite independently and the gap 
remained. Only recently it was realized [TE[] that a class of problems regarding shrink- 
ing domains of Ising spins could be understood in terms of the integrable six-vertex 



model in one of its critical manifolds f[9" |. Besides, using a master equation approach, 
a larger class of nonhermitian and integrable g-deformed models were shown to ap- 
pear naturally as the time evolution operators of several react ion- diffusion processes 



2C 



Within this spirit, our aim in this paper is to find all three-state integrable quan- 
tum chains whose Hamiltonians are time-evolution operators of diffusion-reaction 
processes and whose spectra are known or can be easily calculated. Our motivation 
is twofold: the equivalence of the spectra guarantees us the equivalence of the phase 
diagram and the physical behavior of chemical systems is then in principle deter- 
mined. Consider the long-time behavior of systems governed by the vacuum-driven 
reactions given above, as an example. The systems will relax to a final state where 
there are no particles left, the mean concentration for particle A (or B) decaying as 

CA,B ~ | l-t/r (1-10) 

where a and r are characteristic of each problem. Borrowing the field-theoretic 
jargon we can talk of a massless and massive phase if by these we mean the behavior 
depicted above. In the one- dimensional annihilation and coagulation models one has 
a = \ which imply in a slower decay towards the vacuum state. The inclusion of 
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certain kinds of processes (or reversible reactions) leads to a local steady state, thus 
driving the system towards a quicker (exponential) decay-regime. We shall see in the 
next sections that these behaviors correspond, in the quantum chains to which these 
models are related, to the massless and massive regimes respectively. 

Our second motivation comes from the fact that the connection with integrable 
systems provides us with methods of calculating physical quantities exactly. In this 
respect, we were particularly motivated by the work of Gwa and Spohn |21| who 



calculated the dynamical scaling exponent of the discrete Noisy Burgers Equation. 



This equation was introduced by Burgers |22| as a model for turbulent flow. The 
discretization leads to an asymmetric two-state diffusion problem described by the 

k 

reaction A + # + A where the rates k and r are different. The dynamics of this 

r 

problem is governed by the Hamiltonian of the six-vertex model. 

We summarize the results of this paper in what follows. We found that a large class 
of chemical reactions composed of several simultaneous processes can be understood 
in terms of two quantum chains, namely the L/ g S'[/(2)-invariant model and the 3-state 
f/ (? S'?7(P/M)-invariant Perk-Schultz models. Using the standard basis of matrices 



(E kl ) a ^ = 5k ia Si : b the first chain reads |23[] 
H' = H' + H[ 



L-l 

i+1 



i=l 

+w (e° + e° i+1 ) + a(e° i -e 



L 

-.z 



i=l 

E ll + fizz e z = £11 _ £ 



i=l 

_ mil i 77-22 _z pill 77i22 



;i.ii) 



where Hq is U q SU (2)-invariant and the symmetry-breaking H f Xl which commutes with 
Hqj acts as an external field and does not spoil the integrability of the model. H' Q has 
the same spectrum (apart from degeneracies) as that of the spin-| Heisenberg model 
with an external field and a surface term 

H XXZ = -~ %{at&U-i + °Mn + A <^+i + h + ^ + i) + a (a\ - + (3 

(1.12) 

The phase diagram is well known [24]: without external field, the model is massless 
and conformal invariant when — 1 < A < 1, massive with a ferromagnetic ground 
state for A > 1 and massive with an antiferromagnetic ground state for A < — 1. 
When h ^ the system is massive commensurate for A > 1 — h and massless incom- 
mensurate otherwise, the line A = 1 — h corresponding to a Pokrovski-Talapov (PT) 



phase transition |25! between these two regimes. Since the spectra are equivalent, so 



are the phase diagrams. 
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To see the connection between H' Q and H xxz we rewrite the latter in the basis of 



E matrices 



L-l 



H 



xxz 



8=1 



i+1 



+w (El 1 + EH,) + a {El 1 - E Hi) + b 
where we redefined the parameters of the Heisenberg model as follows 

A + P 



-2A, 



w 



A + h, 



+ h) 



(1.13) 



;i.i4) 



H' is obtained from H xxz by adding extra terms proportional to E^Ef®, and 
Ef^E^ which do not affect the spectrum and correspond, in the chemical scenario, 
to the exchange of the lattice configuration 50 to 05 and 05 to 50 at sites i and 
i + 1 respectively. The diagonal terms are extended accordingly. The wave functions 
of H'q have been calculated but the effect of H[ on the phase diagram has not yet 



been studied ||23|| . We shall reinterpret the known phase-diagram in the language of 
chemical reactions. 

The second class of chains which appear naturally as time-evolution operators of 
chemical systems are the Lf ? S77(P/M)-invariant Perk-Schultz models (PS models). 
The general iV-state Hamiltonian of the PS chains can be written as 



L-l 



H{e 1 , — ,e P+M } 



£ uf' M) 



L-l 

E 



q + q 



-l 



q-q 



2 

-1 N-l 



N-l 



£ EfE% 1 + q + q ' 



-1 N-l 



7+1 



with e = ei 



ep-i 



paa dick 

^ a=0 

£ «^ (a - /?) J^JSf&l } (1.15) 



-€ P 



-ep+M-i = 1. Here (P/M) 



stands for the entire partitions of AT. These models were first introduced by Suther- 



land [2f| for q = 1 and then extended by Perk and Schultz [£7] to q ^ 1. Later |28 
it was noticed that the Uj are generators of the Hecke algebra H(n) (n = L — 1 
defined through 



UiUi 

U t U l±1 Ui - u % 
UiUi + j 



(q + q~ l )Ui i 
Ui±\UiUi±i — Ui±\ 

j > i + 2 



l,2,---,n 



(1.16) 



where to each partition (P/M) of N corresponds an additional set of relations be- 
yond the ones given above, defining the so-called quotients of the algebra. From a 
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mathematical point of view, this underlying algebraic structure has important con- 
sequences as regards the spectrum of each chain j^J. From the physical one, each 
quotient corresponds to systems with different properties. For P + M = 2, one has 
two chains: the (2/0) chain which is the spin-| XX Z model and the (1/1) chain 
which has also been extensively studied [BO] and is believed to describe the coverage 



dependence on fugacity for xenon adsorption on copper [31]. It has a PT phase tran- 



sition in its phase diagram at q + q^ 1 = 2 between a commensurate ordered phase and 
an incommensurate one with oscillating correlation functions. For P + M = 3 the 
(2/1) chain has received most of the attention due to its relevance in understanding 



Anderson's t — J model [32 



We also found a very remarkable property of the mapping of chemical systems 
given by combinations of rates ( |1 . 1|) through (|1.9Q onto the quantum chains we pre- 
sented. It turned out that due to the structure of the mapping not all the parameters 
which we started with (the different reaction rates) are important in settling the 
phase structure of chemical reactions. This will become clear in the next sections, 
were we present our results with more detail. 

The third important point is that for certain chemical models it is necessary 
to define a nonlocal similarity transformation in order to rewrite the Hamiltonian 
describing them in the standard form of the PS Hamiltonians. This kind of transfor- 
mation was already known to map the periodic spin-- XXX with an XY interaction 
term onto the XXZ model with general boundaries proportional to the volume of 
the system [|33||. The transformation we found generalizes it to higher-state chains 



and also takes account of more general Dzialoshinsky-Moriya-type interactions [34| 
parametrized by different variables. 

This paper is organized as follows. In the second section we introduce the formal- 
ism of the master equation on lattices and its mapping onto nearest neighbor quantum 
chains. The strategy we will adopt to identify chemical reactions with quantum chains 
is explained. Section 3 is devoted to the {7 g S77(2)-invariant model, where it will serve 
as a workbench for the application of the ideas developed in the previous section. In 
section 4 we rewrite the PS models with external fields as reaction-diffusion Hamilto- 
nians, constructing explicitly the similarity transformation which map the chemical 
systems onto them. In section 5 we study this transformation and generalize it in two 
directions: first to higher-state models parametrized by different diffusion rates and 
second showing that diffusion-processes on periodic lattices are mapped on chains 
with more general boundary conditions. Finally in section 6 we summon our results 
and some questions we are still faced with. 



2 The Master Equation and Quantum Chains 

The master equation governs the evolution of the probability distribution of 
Markov processes. Due to its almost universal range of applicability, it is one of the 
most important equations of Statistical Mechanics. Here we shall apply it in the 
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context of chemical processes on a chain. Consider a one-dimensional lattice with 
L sites and open boundaries. At each site j we define a variable f3j which takes N 
integer values (0, 1, 2, N — 1). To each possible configuration {(3} = {Pi,---,Pl} 
of the lattice realized at time t we attach a probability distribution P({P},t) whose 
time evolution is given by the master equation 

L-1 



L-l' . 

+ E rS5 +1+m P(/? 1 ,...,A + /,/? fc+1 + m,...,^;t) (2.1) 

l,m=0 ' 

Here and henceforth the prime in a sum over I and m indicates exclusion of the 
pair I — m — 0. The r°'^ are transition rates which equal the probability that, 
in a unit time step and at any site j, a state (Pj,Pj + i) = (a,b) changes to a state 
(P'j,Pj +1 ) = (c,d). We will assume throughout this paper that the transition rates 
depend only on links (nearest-neighbor interaction) and are homogeneous, i.e. site 
independent. The Q a £ are related to the probability that a state (a, b) will not change 
after a unit time step. From conservation of probability we see that they satisfy 

V a , b = (2-2) 

r,s 

With the help of these definitions and the matrices E kl defined in the introduction 
we can rewrite the master equation Q2.1] ) as 



dt 

if we identify \ip) as the probability P({P},t) and H as 

H = £ = £Vi - Ti) 

J=l 3=1 

N-1 

a,b=0 
N-1 ' 



--HW (2.3) 



Tj = Y. r^Bf^: (2.4) 

a,b,c,d=0 

The key to the whole process of identifying chemical processes with known quan- 
tum chains lies in finding the appropriate set of rates so that we can recast H as 

L-1 

H = H + Y,(hi + h i+i +9i~ 9i+i) + Hi (2-5) 
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so that the first two terms on the r.h.s are equivalent to a quantum chain plus external 
field and surface term, and the spectrum of H is independent of Hi. The surface term 
is an extra degree of freedom that we have: since they do not alter the bulk properties 
of the system, we can always define them. We shall look only for those chemical 
processes whose spectrum is equivalent to that of some quantum chain, therefore 
guaranteeing the equality of the phase diagram. However the wave functions are not 
the same. A priori one can define a similarity transformation of the form A(X) such 
that ^ ^ 

A(\)HA-\\) =H Q +Y,{h l + h i+ i + 9i - g i+ i) + H[(X) (2.6) 

such that for some given value A = A we have H[(\ ) = 0. If such a transformation 
were found, then the wave functions of the chemical problem could be obtained from 
that of the quantum chain. No solution has yet been found and the problem remains 
open. 



3 The U q SU(2) model 



With the formalism developed in the last section, we shall now address the 
problem of finding the set of chemical reactions associated to a given quantum chain. 
Our ultimate goal is to use the phase diagram of the chain to explain the chemistry 
of reaction-diffusion processes. In what follows we identify our particles with A = 1, 
B = 2 and inert state = 0. 

We consider a system in which particles A and B diffuse to the right and to the 
left symmetrically, with rates equal to unity 



A 
B 



+ A 
+ B 



rate T 



1,0 
0,1 
^2,0 
0,2 



r 0,l 

1 1,0 
r 0,2 

1 2,0 



(3.1) 



In addition to these processes, we also allow the particles to react according to the 
following 24 vacuum-driven rates: 
• annihilation 



' A + A - 


^0 + 


r 


B + B - 


^0 + 


r 


A + B - 


^0 + 


r 


B + A- 


^0 + 


r 



1.1 

0,0 
2,2 
0,0 

1,2 
0,0 

2,1 
0,0 



(3.2) 
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• coagulation 



death 



polymerisation 



• trapping 



A + A 
A + A 

B + B 

B + B 



A + 
$ + A 

5 + 

+ 5 



A + 0- 


^0 + 


r 


+ A- 


^0 + 


r 


5 + - 


^0 + 


r 


+ 5 - 


^0 + 


r 



A + A - 


^5 + 


r 


A + A- 


^0 + 5 


r 


B + B - 


^ A + 


r 


B + B - 


^0 + 5 


r 



A + B 
B + A 

A + B 

B + A 

A + B 

B + A 

A + B 

B + A 



A + 
+ A 

5 + 

+ 5 

+ A 

A + 

+ 5 

5 + 



r i,i 

1 i,o 
r i,i 

1 0,1 

r 2,2 
1 2,0 

r 2,2 
1 0,2 



1,0 
0,0 
0,1 
0,0 

2,0 
0,0 

0,2 
0,0 



1,1 
2,0 
1,1 
0,2 

2,2 
1,0 

2,2 
0,2 



pl,2 
1 1,0 
F 2,l 
1 0,1 



1.2 
2,0 

2.1 

0,2 

1,2 
0,1 



r 2,l 
1 1,0 



1,2 
0,2 

2,1 
2,0 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



With these processes and the technique developed in section 2 we obtain the following 
Hamiltonian 



H = H + H 1 
Ho = E{-(^m 

i=0 k 



, I r o,i 

^1 1 0,0 



TT.10 pOl 



i)^°^i + (r^ + i)^ 



00 z?22 
i+1 



00 

i+1 
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+(rS;S + r?; + rg + 1% + r$j)j5fJ5, 
+(rio + ril + ^ + r^ + r^)^ 



in 

'i+l 



i22 

'i+l 



ill 

'i+l 



, f r 2,2 , r 2,2 . r 2,2 . r 2,2 . r 2,2\ rn22 tt.22 \ 
+ ( 1 0,0 + 1 0,1 + 1 1,0 + 1 0,2 + 1 2,0 )A ^i+lj 



L-l f 

u J F ' 1 ZT'OO ttiOI , F l,0 ttiOI 77.00 i F 0,2 pOO 77.C 

H l = 1^ \ 1 O.O-^i -^i+l + 1 O.O-^i -^i+l + 1 0,0-^i &i 




_ 1 _p2,l tt.22 77.01 , r 2,l 77.02 77.21 1 r l, 2 77.21 77.02 \ 
+ 1 2,0-^i + 1 0,2-^i + 1 2,0-^i ^i+lj 



(3.7) 



The spectrum of if coincides with that of ifo- We found this solution numerically, 
by successive trials: we started out with a minimum set of reactions, namely those 
corresponding to diffusion and we added vacuum-driven rates one by one until we 
had a 2-body Hamiltonian structure resembling the structure of the quantum chain 
we were interested in. On top of these, we added more processes until we reached the 
point were the spectrum obtained was different. This type of property is characteristic 
of nonhermitian phenomena and has an important physical implication: since the 
spectrum does not depend on each parameter of Hi independently but only on their 
combinations which appear in Hq, our system has 8 effective parameters, which are 
the 8 sums inside parentheses of equation ( |3.7| ) 2 . Yet, the problem of identifying 
this effective chemical Hamiltonian with the chain given by H' in ( |1 . 1 1| ) requires the 
reduction to an even smaller set of parameters since H' depends only on four. The 
solution is to find the proper way of recombining the eight effective parameters of the 
chemical model into a final set of four. We will solve the problem in steps. 

First, to avoid carrying too heavy a notation throughout the text we rename 
groups of rates regarding the same category of processes as follows (i = 1,2) 



2 Due to conservation of probability each diagonal element of a Hamiltonian must equal the sum 
of the nondiagonal entries in that same column. We can see this if we look at the definitions of rates 
in section 2 and realize that the conservation of the probabilities implies in a relation among rates 
which can be written as X^fc^z Hki = Hu. 




.2,1 
0,0 



annihilation 



annihilation 



9 



Df = T ' l ± T^o death 

Cf = T*' + Tq \ coagulation 

Pi = r ' 2 + rl 'J polymerisation 

P 2 + = Tg'^ + T^'o polymerisation 

= [rlj ± r$i] + [rj;? ± rJJ] 

+ [ r io ± ^2] + [ r o;2 ± rjo] trapping (3.8) 

The problem is solved through the rearrangement of the rates into the following 
variables 

v = A x + C+ + P x + - D+ - 2 

„ = ^±^ + 1 

y 4 

a = T~ 

l = (3.9) 
together with 3 conditions on the rates 

2w + 2g + v = A 2 + C 2 + + P 2 + 
Aw + 2v = A+ 2 + T + 

Al = A u +T- (3.10) 

After some simple algebraic manipulation, we can rewrite Hq as a new Hq which 
reads 

Ho = £ " + El°E^ + EfE^ + EfE^) 

i=i 

+w(e° + e° i+1 ) + g{e\ + e z i+1 ) + a(e° - e° i+1 ) 
+l(e z -e z i+1 )+ve° i e° i+1 
e° = E 11 + E 22 e z = E 11 — E 22 (3.11) 

Comparing this expression with of eq. ( |1.11|) we conclude that they have the same 
phase diagram since the term in H having the parameter / as coefficient is a surface 
contribution. To interpret the phase diagram of the chemical model in terms of the 
phase diagram of the XX Z model we still have to recover the ?7 g S , ?7(2)-symmetric 
spectrum, which can be done by requiring that g = I = 0. In terms of chemical rates 
this condition means 

Dt = E+ 

£>r = D 2 (3.12) 
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This amounts to saying that A and B are indistinguishable as we can see by examining 
equations (|3.9|), ( |3.10| ) and ( |3.12| ). It suffices now to identify the parameters of the 



Heisenberg chain with the rates of our chemical model. We obtain 

A 1 + C? + P+ 



2 

Df - (A l + C+ + P+ 



h 



A = 1 

2 

Dx 

a = 

P = -!-^ + C ' ++ /' + + C ' (3.13) 

The analysis is straightforward. Making h = implies that no rates but death 
survive. In this situation 

2 

13 = -A (3.14) 

By varying the death rate we go from a massive regime to a massless one, that is the 
time evolution for the concentration of particles has an exponential or algebraic fall- 
off respectively. This can be understood on physical grounds: being a 'one-particle' 
process, death happens irrespectively of any other processes ocurring in the system, 
i.e. it is not diffusion-limited since any particle can die alone. It therefore outruns 
the characteristic time scale set by diffusion and brings about a quicker decay. With 
an external field we have 

n+ 

A + h-l = ^- (3.15) 

If we take death with a non-zero probability, then from the equation above one sees 
that our system is massive ferromagnetic. In Chemistry this means having a ground 
state with no particles (we identify spin up with 0) By varying the rate of death we 
approach the line on the phase diagram given by A + h = 1 . The system undergoes a 



PT transition |25| when the energy of the state with just one particle equals that of 
the state with no particles and it becomes the ground state. We have a level-crossing: 
since death is absent the system can evolve to a final steady state where only one 
particle is left. 

4 The Perk Schultz chains 

4.1 The [/ g S^(P/M)-invariant chain 

It is a well established fact that the Hamiltonian of the 6- vertex model is the time- 



evolution operator for the two-state asymmetric diffusion process [18|. On the same 
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ground we expect that the PS chains will play the role of time-evolution operators of 
higher-state asymmetric diffusion processes since they are the Hamiltonians of higher- 
state ice models. This picture is however far from complete and can be extended to 
encompass more general chemical systems if we reinterpret the additional reactions 
on chemical Hamiltonians as external fields in the PS chains they are mapped onto. 
For the sake of completeness, we present first the results without external fields and 
then we proceed with the more general models. 

Following we first consider a system in which particles A and 5 diffuse and 
interchange positions on the lattice according to 



+ A 
5 + 
+ 5 
B + A 
{ A + B 



A + 
+ 5 
5 + 
A 
B 



B 
A 



rate T r 



(4.1) 



With these processes and the rates defined above we get a Hamiltonian which reads 



H 



.7=1 k 



11 + E m E 22 



7-.ll T-.22 N 

h 3 h S+l< 



+V R {EfEf +1 + EfEf +1 + EfEf +1 ) 
-V r ( E 10 E m + E 20 E 02 + E 21 E 12 ) 

-T R {EfEf +l + EfEf +l + EfEf +l )) 



(4.2) 



Defining \ff^ = q, which measures the asymmetry of the diffusion, and the diffusion 

constant \/T L T R = V, which sets the time scale of the problem, this can be rewritten 

as 



H. 



V 



+ 



q + q 



I' 1 E E f E Hi + 9 E «i 



a>/3 



a</3 



-1 2 



X; e a Ef a E^ + 



q-q 



-1 2 



+1 



a=0 



(4.3) 



where (eo,ei,e 2 ) = (1,1,1). Comparing with eq. (|1.15|) we see that this is the 
U q SU(3/0) chain in a non-standard form. To bring it to the standard form we have 
to consider a nonlocal similarity transformation S which is given by 



ai,a 2 ,--,a L =0 



E 



OtL<XL 



(4.4) 



The implications of this transformation coming from boundary terms are discussed 
in detail in the next section. If we apply it to the Hamiltonian (|4.3| ) we obtain the 
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standard U g SU (3/0) chain which reads 



+ e ^r^ a i + e »w« - mr E llx 

Z a=0 Z 



4.5) 



We use two notations to differentiate between what we call the non-standard (non- 
hermitian) and standard (hermitian) forms of the chain. 

To reproduce the {7 9 S77(2/l)-invariant PS model, two new processes have to be 
added, namely those corresponding to coagulation of B 



B + B -> B + rate T 2 2 ' 2 Q 

b + b^® + b rjg ( 4 - 6 ) 

such that r^o + Tq^ = T>(q + q^ 1 ). The main difference now in comparison to the 
first case we studied is that we will have an extra diagonal contribution equal to 
r^'o + Tq' 2 an d two new nondiagonal pieces proportional to each of these two rates. 
The Hamiltonian that we find in this case reads 

it TT L—l f p2,2 -p2,2 

n _ -"{1,1 -1} J i 2,0 ^22 pi02 1 0,2 ,-,22 pi20 



V V 



t{^EYEf +l -^EfEf +1 ) (4.7) 



where H^^-iy is obtained from ( |4.3| ) with (eo, 61,62) = (1,1,-1). The spectrum of 
H is equal to the spectrum of if^i-i} and again we can study only the properties of 
the latter as long as we are interested only in the phase diagram. H^i-iy is the non- 
standard representation of the (2/1) PS chain. With the similarity transformation of 
eq. (O) we obtain the standard L/gS77 (2/1) -invariant PS model 



SH^iyS- 1 = # (2/1) (4.8) 

where if^ 2 / 1 ) is given by (|1.15| ) with (eo, 6i, 62) = (1, 1, —1). 

Finally, the last PS chain can be reproduced by the inclusion of the rates corre- 
sponding to coagulation of A 

A + A -> A + rate r\'l 

a + a -> + a r$ ( 4 - 9 ) 

with T{q + = P(g + g _1 ). The Hamiltonian in this case reads 

TT TT L— 1 , pl : l plil n 

-" _ -"{1,-1-1} J 1 1,0 pll F 01 1 0,1 pill plO 1 (AX(\\ 

v - ~~d — ~ ~w E i (4 - 10) 
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The same argument applies: the spectrum of H is equal to the spectrum of Hjx,-i,-i}, 
which is the non-standard L r ? S'L r (l/2)-invariant PS model. This is in turn equivalent 
to the standard form through the trasformation S 

SH^^S- 1 = # (1/2) (4.11) 

We would like to point out one particular feature of the two last chains. Once we 
find the spectrum of the (2/1) chain, the spectrum of the (1/2) can be obtained in a 
straighforward manner: we reverse the sign of the whole spectrum of the (2/1) model 
and an overall constant equal to (L — l)(q + q^ 1 ), which is the highest energy of the 
(2/1) chain (this value is indeed the highest energy for all L-site (P/M) chains with 
q real and non-zero P and M. An outline proof is given in ref. ||29|| ). In the chemical 
scenario however, they correspond to different physical pictures and the reason is the 
positivity of the spectrum. 



4.2 The U q SU (P/M)- invariant chain with external fields 

We now generalize the results above to include the most general set of vacuum- 
driven processes which can be written as a pure nonhermitian Hamiltonian which 
does not contribute to the spectrum of the PS chains. We analysed this problem on 
the computer and found that beyond the 24 rates of section 3 we can include one 
more set of reactions, namely those corresponding to mutation and transmutation of 
A. They are defined through 



A + 0- 


-> £ + 


rate V 




^0 + 5 


r 


A + 0- 


-> + B 


r 


+ A- 


^5 + 


r 



1.0 
2.0 
0.1 
0.2 

1.0 
0.2 

0.1 
2.0 



(4.12) 



For these reactions we define the following combination of rates 



r o,i , r i,o 

1 0,2 ^ 1 2,0 

r o,i , r i,o 

1 2,0 ^ 1 0,2 



(4.13) 



The problem is the same of section 3: the original chemical process we start with has 
more parameters than the PS chains. The question is whether it is possible to reduce 
these rates to the appropriate set of parameters in a way that these can be rewritten 
as external fields and surface terms. We found a positive answer in all three cases 
studied. Using the master equation approach we get a Hamiltonian which can be 
written as two separate pieces as follows 



H 
V 



H 



{f*} 
V 



L-l 



]T (hi + h i+1 + gi - g i+1 ) + 



(4.14) 
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where H{ Ea y is given by (|4.3| ) and H^> = ^j-ffj being the same in all cases, reads 



H, 



(i) 



+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 



r 0,l rpOO 
p2,0 T-.02 



77.00 , 



e: 



01 
i+1 
22 



2,2 77.12 77.02 



IW 2 



3+1 + 
^+1 + 



1.2 77.01 77.12 



1 0,1-^ 

p2,l 77.22 

1 2,0-^i 

pl,0 77.21 

1 2,0-^i 



pl,0 77.01 77.00 , r 

1 0,0-^i + 1 

pl,l 77.01 77.01 , p2 

pl,l 77.01 77.11 | p2 
pl,l 77.21 77.01 



F 2,2 7-.02 77.12 

pl,2 77.11 7702 

-1,0^, 

p2,l 77.12 77.01 

-1,0^, Ei+x 

p2,l 77.02 77.21 
pO,l 77.00 77.21 



+ r 
+ r 

+ rj> 
+ r 
+ rj 
+ r n 



2 77.00 77.02 

(A ^i+i 

2 77.02 77.02 

O^i ^i+1 

2 77.22 77.02 

1 77)01 7712 

2&i E i+1 

2 77.01 77-02 
1 77.02 7711 



2 
2 

2 

0"i 

0^-01^-20 



7^1 

j^22j^02 



i+1 



(4.15) 



The /i's and g's are the external fields and surface terms. Since they are different for 
each chain, we will present them separately. Due to the symmetry properties of the 
PS chains, the field (and surface contribution) are decomposed into 3 independent 
variables h a (g a ), one for each conserved particle of type a. Since one of our particles 
is an inert state, we take the corresponding variables h (g ) equal to zero. 

i) The U q SU(3/0) chain 



For the first of our Hamiltonians the field and surface terms are given by 

M+ + Xt + Df 



2V 



2V 



Dl 
2V 

2V 



-9i 



~92 



(4.16) 



together with the conditions 



M+ + X+ + D+ = C+ + P+ + A 1 



Af 2 + T + 
A u + T- 



P 2 + + A 2 
+ 



M+ + X? + Dl - 

D 2 --(Mr + xr 



Dt 



(4.17) 



ii) The U q SU{2/\) chain 
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In this case we get the following relations between field, surface contributions and 
reaction rates 



M+ + Xt + Dj 
2V 

Mf + Xf + Dj~ 
2V 

Dt 
IV 

EL 

2V 



= h 

= -9i 

= -h 2 

= -92 (4.18) 



For these identifications to hold we have to impose the extra set of relations among 
rates 

M+ + T+ + D+ = C+ + P+ + A! 

D+ = Ct + P^ + Az-Diq + q- 1 ) 
A++T+ = M+ + X+ + D+ + D+ 

A~ 2 + T- = D--(M-+X- + D-) (4.19) 

iii) The U q SU{\/2) chain 

For the last model we have 

M+ + Xt + Dt 



2V 

Mf + X£ + Dj 
2V 

Dt 
2V 

£i 

2V 



-9i 
-hi 

-92 (4.20) 



The consistency conditions for these equations are 

M+ + T+ + D+ = Ct + Pt + M -Diq + q- 1 ) 
Dt = Ct + Pt + Ai-Diq + q- 1 ) 
A+ 2 + T + = Mt + X+ + Dt + Dt 

A u +T- = D 2 - (Mf + Xr + £>D (4-21) 
We observed that if in place of mutation and transmutation of A we considered the 
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corresponding processes for the particle B 



' + B -> + A rate V\ 
£ + ^ + A T; 



,0,2 



< 



5 + 0^A + 



+ 5^ A + 



r 



r 1 



,0.2 

1,0 

,2.0 

1,0 



(4.22) 



from which we define the relation of rates 



m ± — r ' 2 4- r 2 ' 1 



"± _ p0,2 i F 2,0 
2 — 1 1,0 31 1 0,1 



(4.23) 



the spectrum also remained invariant, but not if we took mutations and transmuta- 
tions for both particles at the same time. The reason is that the inclusion of both 
processes in the system yields a local steady regime given by the reversible process 
A + # B + 0. In the presence of reversible reactions we expect a totally different 
physical picture which translates itself, among other things, into a different spectrum. 
The phase diagram of the chemical processes can now be in principle explained in 
terms of the Physics of the PS models. Unfortunaltely only the phase diagram in the 
absence of fields is known: this is however not so interesting because it implies then 
that many processes do not survive - and the whole mapping only makes sense if one 
is able to study non-trivial cases. We will return to this point in a future publication. 

5 Similarity Transformation and Boundary Con- 
ditions 

We shall now consider more thoroughly the transformation given by the matrix 
S which we used to rewrite our 3-state Hamiltonians describing diffusion processes 
as the standard PS chains. We shall address here two points: the first regards the 
extention of S not only to higher-state models but also to multi-parameter diffusion 
processes, i.e. those characterized by a larger set {q} = {q±, q 2 , • • •} of parameters in 
substitution to the one-parameter diffusion we dealt with so far. Second, we want to 
look at the effect of S on periodic chains and the physical features it induces. 

5.1 Multi-parameter and Higher-State Diffusion Processes 

As before we consider a system in which (N — 1) particles diffuse with and 
interchange places among themselves. We define the following rates 



A x + A. 



-y 



Ay + A x r 




(5.1) 
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from which we also define N( - N 2 ^ parameters q xy 



\ y a ^ a I = qxy ^' 2 ^ 

\ Ay, A X 



with q xy = q yx x . These are the generalization of the parameter q defined in the 
previous sections. We impose a homogeneous time scale 



/ r 0,Ai r Ai,0 _ / r 0,A 2r A 2 ,0 _ _ r A N _2,A N _ 1 - r A N _ 1 ,A N _ 2 _ ^ /r o\ 

V 1 Aj.O 1 0,Ai ~ V 1 A 2 ,0 L 0,A2 - ' " _ V ^iv-i,A JV -2 i A n _ 2 ,An-i ~ U 



and obtain a nonhermitian Hamiltonian which reads 

jj L-l JV-1 
k=l ^Xy^y=0 



TT L-l ( N-l s 

p = E { E &„ (>r^i - ^r^O } (5.4) 



We want to find a similarity transformation which makes H hermitian. The (nonlo- 
cal!) similarity transformation which accomplishes this reads 



N-l , N-l s 

S = Y. ( II (q X y) fxy{ai,a2 '"' ,aL) ) E aiai <g> E a2a2 (8) ... (8) (5.5) 

ai,a2, •••,(*£ =0 v x>j/=0 

with the functions given by 

1 L 

fxv(ai,a 2 ,---,oiL) = g E (^n,^^ - 6a n ,Am,*) (5.6) 

n>m=l 

We point out in the expression for the function that the order of its arguments is 
very important since they index the sites on the chain, that is /(•••, otj, • • •) ^ 

/(•••, ttj+l, Q!j, •••). 

Proof 

For clarity, we will consider the diagonal and nondiagonal pieces of the Hamilto- 
nian separately. Since S(J2j ^jj+i)^ -1 — Z)j SHjj+iS -1 , it suffices to consider the 
action of 5 on the 2-site operator only. 

Consider first the action of S on the diagonal piece of the Hamiltonian. We have 

j^diag N-l , 

S S" 1 = II ( q^ ai ' a2 ''"' aL) E aiai <g> E a2a2 (8) ... (8) J 

N-l 



x ( E Q*E?E» +1 ) 

N-l , \ 

X E II (q-^^'^-'^E 0101 (8) (5.7) 



{/3}=0 
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We rearranje the terms in the expression above as follows 

Tr(diag) N—l N-l N-l N-l , 

S flM±Ls-' = e e n n (g^^"^ q -t^^^ 

W {(}} x>y=0w>z=0^ 

x E aiai E^ 1 ® E a2a2 E^ 2 ® ■ ■ ■ ® <g> <g> ■ ■ ■ ® E aLaL E pLpL 

N-l , x 
x E E [qabE ajClj E aa E^^ ® E^+ ia ^+ 1 E bb E^+ lP ^ 1 ) (5.8) 

Since re, u> and z serve only to index the same set {q} of diffusion parameters, we 
can take x = w, y = z. Also, the multiplication properties of the matrices E a/3 gives 
us 

E pq E rs = 5 q:r E ps (5.9) 

We have therefore 

Tj-(diag) jv-l N-l N-l , 

S S' 1 = E E II ( q f x xv{ai ' a2 ''"' aL) ~ fxvi ~ l3lA ''"' (3L) 

V {a} {(3} x>y=0^ 



x E aiai 5 aiA ® E a2a2 5 a2A ®---® ® ®---® E aLaL 5 aLt(3L 
x E E (qabE] a E bb +1 5 aj>a 5 aj+ub 5 aj>(3j 5 aj+1 ^ +1 ) (5.10) 

.-.a.- i 1 „-/-h— n \ ' 



Taking into account the Kronecker's deltas we have 

jj(diag) N _i A r_ 1 A r_ 1 

s ^±L s -i = E II qiir {ai ' a ^-' aL) - f ^ ai ^^ aL) [ E qabE* a Ef +l \ 

{a} x>y=0 ^a^b=0 ' 

x E aiai ® E a2a2 ® ■ ■ ■ ® ® l (i+1) ® ■ ■ ■ ® E aLaL (5.11) 
The exponent of q xy is clearly zero. Summing over {a} we have finally 



S^^S- 1 = E 5a*l®l®-"®£$®£g+l) 



(5.12) 



This concludes the proof that the diagonal elements of the Hamiltonian are not 
changed by the transformation generated by S. 

To see the effect of S on the nondiagonal piece we proceed as before. We consider 
only the 2-body operator acting on sites (j, j + 1) and we especialize to one given 
value of the pair (a, b). We have 

S^L^S-i = s(q ab EfEf +1 + q^EfEf^^S- 1 (5.13) 
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As in the diagonal case after multiplying the E a/3 matrices on each site we obtain 

Tj(nond) jv_! JV-1 , N 

^ {a} x>j/=(A ' 

x E aiai 5 aiA ® • • • (8) (8) <g> • • • <g> E aLaL S aL , Pt 
x ^ 



XI ( QabEfEf +l 8 a . ta 8p. >b + qjEfEf +1 5 a] fi 5 p . >a 



a>b 

X tfoH.^A^+i ( 5 - 14 ) 



This reduces to 



Tj{nond) jv-1 JV-1 

g g-1 _ J~J ^/ sv (o!i,-,aj,a^+i,—,a£)-/ < E»(aii,-,aj + i,aj,—,az,) 

{a} x> 2 /=0 

X <g> E a2a2 ® ■ ■ ■ <g> <g> ® • • • ® J7J a LOL 

x i .. uabipab x X i „-l IT'&a 17x16 



q ab EfEf +l 5 a3A 5 a3+l , b + q- ab l EfEf +l 5^ b 5 aj+lA (5.15) 

We have to analyse the exponent of g xy now. We once again point out that the order 
of the arguments of the functions / is important. With this in mind we have 

f X y(oci, ■ ■ ■ ,aj,a j+ i, ■ ■ ■ ,a L ) - f xy {oii, ■ ■ • , aij+i, atj, • ■ ■ , a L ) = 

^Oj+i ,x&a.j ■ ,y ^c«w_|_i ,y^oij ,x (5.16) 

Substituting this expression on (5.15p we get 



rrinond) jV-1 A r_ 1 . 



^ {a} x>y=0 

x <g> jT} a 2 a 2 g . . . g ]_(?) g g . . . g (5.17) 

Summing over {a} gives us an identity matrix on each site. The deltas also 'kill' the 
sum over aij, atj+i such that we will be left with 

g ff' 1 = [] [q S ^ Sa > v ~ Sa > xSb ' v qabEfE]l x + q^ 

U x>y=0^ ' ' a > h 



Ej E j+1 + q ab q ab E j E j+1 
{EfEf +l + EfEf +1 ) a>b (5.18) 
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This concludes our proof. 

To finish we would like to indicate how in the one-parameter iV-state diffusion 
process this expression simplifies. The rates are defined according to 

A x + A y ^A y + A x T Rl T L x>y = 0,l,---,N-l (5.19) 

for which the L-site Hamiltonian reads 



' ' V Z L a>0=O a<0=O 



k=l 



(5.20) 



This is the non-standard ?7 (? S'f/(A^/0)-invariant PS model. To derive the standard 
form we consider the matrix S from eq. ( |5.5|) with q xy = q for all x, y. By observing 
that in this case we have 

1 L 

fi(a 1 ,a 2 ,---,a L )-\ h /zv-i(ai, "2, • • • , ot L ) = - sign(a n -a m ) (5.21) 

n>m=l 

we get for S the following expression 

N ~ l 1 L 

S = ^ En>m=i «<ffn(«n-a») g, £a 2 a 2 g . . . g ( 52 2) 

It is a simple algebraic exercise to bring the chain ( |5.20| ) to its standard form under 
the action of S. Actually the result holds for any of the f/gS , f/(P/M)-invariant PS 
chains with P + M = N. 



5.2 Periodic Boundary Conditions 

Due to its nonlocality, the effect of S on free or periodic chains is different. In the 
first case it brings the original non-hermitian chemical model to a hermitian quantum 
chain. In the second case, it brings the non-hermitian chemical model to a hermi- 
tian quantum chain but with general boundary terms which imply in a generalized 



Dzialoshinsky-Moriya-type of interaction in the bulk [p4}| . This transformation was 
known for many years for 2-state systems P3| . 

To understand the effect of S on periodic chains, it is appropriate to start with 
the simplest system, namely a particle A diffusing to the left and right with rates 

A + dS^dl + A T Rl T L (5.23) 

where T> and q are defined in the usual way (see section 4). The Hamiltonian is given 
by 

P = E{ — 2 — ^ k k+1 2 — ^ stgn{a-(3)E k E p k p +1 

k=l a=0 ~ a^0=O 

-^E^E^-qE^EllA (5.24) 
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It is convenient to change to the more familiar basis of Pauli matrices, by making the 
identification 



E 00 

E 11 

E 01 
E 10 



l + a z 
1-a 2 



a 

a' 



(5.25) 



In this basis the Hamiltonian can be written as 



H 



q-q- 



(5.26) 



The constant term is irrelevant in our discussion and we will drop it for the time 



being. Applying eq. (|5.22j ) to the remaining terms in the chain we get 



s-s-i = * 

V V 



Z k=l 1 Z 
-\{2q L aia^ + 2q- L a L a+ + ^f^^l 



(5.27) 



We can bring this under a common summation sign if we define the boundaries as 



'L+l 



a 



L+l 



a 



L+l 



q L °t 



(5.28) 



With this definition we have 



H lJ^r _ q + q- 

P = 2^ \ 2(T k °k+i + 2o k <7k+i + — 2~ 



-i 



-o h o 



k u k+l 



(5.29) 



Since q is real, we have a boundary term which is proportional to the volume of the 
system, changing the whole structure of the problem. In other words, the dynamics 
of diffusion processes in a chain with open boundaries are different from that of a 
periodic chain |35| . 

For higher-state diffusion models on periodic lattices we will get, after applying 
S, chains which cannot be put in terms of simple generalized boundaries. Rather 
we get a bulk interaction. To see this we first note that the whole effect of S is to 
change the off-diagonal terms of the 2-body chemical Hamiltonian at the boundary, 
namely Ef E\ a and E b L a Ef. In the 2-state case there are only two possibilities for 
the pair (a,b): (1,0) and (0,1). Since opposite pairs are related by an inversion of 
the power of q, we have effectively one result which is q L . In higher-state cases, the 



22 



power depends on the particular (a, b) chosen. We present the results for the 3-state 
model and the for the general TV-state case in what follows. 

i) 3-state diffusion process 

Here we have (a, b) = (0, 1); (0, 2); (1, 2) and the reversed pairs. The action of S 
yields 



^{g- 1 ^ 01 ^ }^- 1 = q~' 2 Ef <g> Mi <8> M x <g> ■ ■ ■ <g> M 1 ® E% 



-2 7T.01 



10 



Mi = 



V 



1/ 



(5.30) 



for (a, 6) = (0, 1). Exchanging a, 6 we get the same structure but with q replaced by 
q^ 1 as expected. For the pair (a, b) = (1, 2) we obtain 



Siq^E^Ef }S~ 1 = q- 2 E{ 2 ®M 2 ®M 2 

(I \ 



M 2 ® 



Mo 



V 



(5.31) 



9/ 



Finally for (a, 6) = (0, 2) we have 

Siq^EfEf }S~ 1 = q- 2 Ef ® M X M 2 ® M X M 2 



M X M 2 ® E L 
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(5.32) 



We see now that in the case under study, the effect of the transformation cannot be 
put in the boundary terms in contradistinction to the 2-state problem - the interaction 
is 'spread' over the bulk. The structure is analogous in higher-state models as it is 
shown next. 



ii) iV-state diffusion process 

Assume that (a, b) = such that i < j — 0, 1, 
S{q- l E^E ] *}S- 1 = q~ 2 El j ® M ® M (g 



case 



• ,N — 1. We obtain in this 



M ® E 3 l (5.33) 
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where M is an N x N matrix which reads 



/I 



M 



(5.34) 



V 



1/ 



where q appears at the i-th and j-th positions in the diagonal. All other are 
equal to 1 for k < i and k > j and equal to q 2 in the case i < k < j. The interaction 
is spread along the bulk as we can see from eq. ( |5.33|) . 



6 Conclusions 

We must start this summary with one question: what is the purpose of this whole 
technique? Does it bring any new information to what is already known regarding 
reaction-diffusion processes? 

Reaction-diffusion models are per se a very rich and fascinating field, from a physi- 
cal and mathematical point of view. The range of physical phenomena which they en- 
compass is extensive, from dispersive transport in amorphous silicon [^] to polymeri- 
sation processes like thermal soliton-antisoliton interaction of trans-polycetylene |p~3| , 



to name a few. Experimental studies have also been conducted on one- dimensional 
coagulation models ]57| . This alone justifies their study and any new method which 
might shed some light into the structure of these problems is highly welcome. 

We showed in this paper that a large class of chemical processes can be understood 
in terms of the U q SU(P/ M)-mvariaiit Perk-Schultz models and the U q SU (2)-invariant 
model with external fields. Since both are integrable, this opens a new perspective, 
namely of employing analytical methods to calculate physical quantities of interest 
exactly. We showed how the phase diagram in the latter case completely settles the 
phase diagram in the chemical model, the corresponding massless and massive regimes 
yielding two types of time-dependence for the decay of the particles' concentration. 
We showed also what processes correspond to a field in the language of quantum 
chains and that not each and every processes governs the Physics but combinations 
of them. 

In the Perk-Schultz models little is known about the phase diagram with external 
fields. Since the effect of a field on the spectrum is trivial, the problem is a feasible 
one. However the limitation comes mostly from hindrances in numerical techniques. 
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To calculate the phase diagram one needs data from large lattice sizes in order to 
make good extrapolations. Here the integrability of the model comes in handy. The 
Bethe Ansatz Equations which allow for good spectral evaluations are not known 
in the case of open chains. For periodic chains we must take care. The chemical 
Hamiltonians associated to the Perk-Schultz models are usually not in a standard 
form. We found a similarity transformation which brings them to the standard form, 
and whose properties were not yet discussed in the literature except for the 2-state 
processes. This transformation induces terms of physical relevance in the bulk which 
are generalizations of the Dzialoshinski-Moriya interaction, thus confirming the fact 
that the Physics of nonequilibrium problems depends on the boundary conditions 
imposed |55| . The three-state Perk-Schultz models with these new interactions are 
soluble through the Bethe Ansatz technique, which in this case yields the following 
set of coupled nonlinear equations 



^(ATo+at^M! ( sinhjX^ +e 7/2) y 

7 /2)/ 



j(Ni +N 2 ) ,M 2 ,M 
fc 2 fc 



n 

a=l 
Mi 

n 

a=l 
M„ 

n 

a=l 
Mi 

n 

Q=l 



sinh(X^ 


_A<°) 

a 


+ ei7) 


sinh(\^ 


x(0) 
— Aa 


- eo7) 


sinh(X^ 




- 6i7/2) 


sinh(X^ 


A (1) 


+ ei7/2) 




- \ (1 > 


+ e l7 /2) 


sinh(Xa^ 


- A,'" 


- e l7 /2) 


sinh(\£> 


- A, (1 » 


- £2 7 ) 


sinh(\a^ 


- A,'" 


+ eii) 



(6.1) 



where {e a } are the parameters of the Perk-Schultz chains, and 7 is related to q 
through q = expiry). The N^s correspond to the particular number of particles of 
type i in each charge sector and the Mj's equal the number of roots of the set of 
coupled equations. They are obtained from the Ni's through 

M i = L-(N + N 1 + --- + N i ) (6.2) 

where L is the lattice size. We shall analyse this problem in a future publication. 
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